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Within a perturbative approach to Quantum Chromodynamics (QCD), we show how 
to extend ordinary DGLAP longitudinal evolution equations to include the radiative trans- 
verse momentum generated in the collinear branching regime. Considering Semi-inclusive 
Deep Inelastic Scattering as a reference process, we perform such a generalization both in 
the current and in the target fragmentation region. These distributions are then used to 
predict semi- inclusive Deep Inelastic Scattering cross-sections onto the whole phase space 
of the detected hadron. 
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1 Introduction 

Semi-inclusive Deep Inelastic Scattering (SIDIS) cross-sections, at variance with 
inclusive DIS, are differential on final state hadron variables and thus are a rich 
source of information on perturbative and non-perturbative QCD dynamics. In 
presence of a hard scale, set by the exchanged virtual boson in DIS, perturbative 
QCD predicts, thank to factorization theorem pQ, the SIDIS cross-sections. The lat- 
ter theorem garantees that non-perturbative distributions which parametrize soft 
hadronic wave functions are process independent and obey renormalization group 
equations, i.e. DGLAP evolution equations [2]- Furthemore final state hadrons in 
SIDIS are expected to have a sizeable transverse momentum generated by QCD 
radiation and this effect is indeed observed in experimental data, see for istance [3] . 
In the following we show how to extend renormalization group equations to include 
transverse degree of freedom. We perform such an extension for space and time-like 
evolution equations in the current framentation 0] E] , where factorization theorem 

strictly applies. In the target fragmentation region a perturbative description 
of the SIDIS process can be achivied through fracture functions 6 , which obey an 
inhomogeneous evolution equations and for which collinear and soft factorization 
has been proven [7J|H]. We thus extend fracture functions to include transverse de- 
gree of freedom and generalize their evolution equations. The SIDIS cross-section 
is then presented at leading logarithmic accuracy in term of transverse momen- 
tum dependent (TMD) distributions. Numerical methods for solving the evolution 
equations and physical issues related to the obtained solutions are discussed. 



Czechoslovak Journal of Physics, Vol. 56 (2006), Suppl. ? 



?1 



2 Trasnverse momentum dependent evolution equations 

Ordinary DGLAP evolution equations |2] resum leading logarithm (LL) of the 
type [(a s /27r) log(Q 2 //i^)]" originating from quasi-collinear partons emission con- 
figurations, where u. 2 F represents the factorization scale. Leading contributions are 
obtained when the virtualities of the partons in the ladder are strongly ordered, 
Q 2 3> fc 2 j_ 3> ■ ■ • 3> k 2 j_ . At each branching, the emitting parton thus acquires 
a transverse momentum relative to its initial direction. The radiative generated 
transverse momentum can be taken into account through transverse momentum 
dependent evolution equations 



One-particle distributions 2?^(z^, Q 2 ,p±) in eq.Q give the probability to find, at a 
given scale Q 2 , a hadron h with longitudinal momentum fraction Zh and transverse 
momentum relative to the parent parton i. Pij(u) are the time-like splitting 
functions which, at least at LL accuracy, can be interpreted as the probabilities 
to find a parton of type i inside a parton of type j and are expressed as a power 
series of the strong running coupling, Pij(u) — ^ n=0 &™(Q 2 )Pij l \u). The order n 
of the expansion of the splitting function matrix Py (u) actually sets the accuracy 
of the evolution equations. The radiative transverse momentum q± at each branch- 
ing satisfies the invariant mass constraint q\ = u (1 — u) Q 2 while the transverse 
arguments of V^(zh, Q 2 , p±) are derived taking into account the the Lorentz boost 
of transverse momenta. The unintegrated distributions fulfil the normalization: 



This property garantees that we can recover ordinary integrated distributions from 
unintegrated ones. The opposite statement is not valid since eq.Q contains new 
physical information. In analogy to the time- like case we consider now a initial state 
parton p in a incoming proton P which undergoes a hard collision, the reference 
frame being aligned along the incoming proton axis. We thus generalize eq.Q to 
the space-like case |1U|: 



One-particle distributions T\,{xb-, Q 2 , &_!_) m eq.(J3J) give the probability to find, at a 
given scale Q 2 , a parton i with longitudinal momentum fraction xb and transverse 




(1) 
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momentum k± relative to the parent hadron. The unintegrated distributions fulfil 
a condition analogous to the one in eq.J5J), i.e. : 



J d 2 k ± r P (x B ,Q 2 1 k ± ) = r P (x Bl Q 2 ). 



(4) 



We note that the inclusion of transverse momentum does not affect longitudinal 
degrees of freedom since partons always degrade their fractional momenta in the 
perturbative branching processes. 

The perturbative approach can also be extended in the target fragmentation 
region introducing new non-perturbative distributions i.e. fracture functions 
Ai Ph (x, z, Q 2 ), whose collinear and soft factorization has been proven respectively 
in [f| |Hj . These functions express the conditional probability of finding, at a scale 
Q 2 , a parton i with momentum fraction a; in a proton P while a hadron h with mo- 
mentum fraction z is detected. Fracture functions are shown to obey inhomogeneous 
evolution equations |SJ . We generalize these distributions to contain also transverse 
degrees of freedom. The fracture functions A4 P h (x, k±, z, p ± , Q 2 ) give the condi- 
tional probability to find in a proton P, at a scale Q 2 , a parton with momentum 
fraction x and transverse momentum k± while a hadron h, with momentum frac- 
tion z and transverse momentum p ± , is detected. Under these assumptions the 
following evolution equations can thus be derived |10|: 

dMp^x, k ± , z, p ± , Q 2 ) _ a s {Q 2 ) j f 1 du D . . f d 2 q 



dQ 2 2tt 



1-z 



V? J 7T 



■mU*, z, P± , Q 2 ) + I* -jJ*WW**L S( (1 - u)Q 2 - q \) 

• \U U ' Jx X(l — U)U Z J TT 

, .. o k i — o i \ _>, / zu 



The homogeneous term has a pure non-perturbative nature since involves the frag- 
mentation of the proton remnants into the hadron h, while the active parton i 
evolve towards the hard boson vertex. The inhomogeneous one takes into account 
the production of the hadron h from a time-like cascade of parton j. In this latter 
term Pjj' (u) stand for the unsubtracted splitting functions 11 . Fracture functions 
fulfil the normalization condition: 

J d 2 k ± J d 2 P± M i Pyh (x,k ± ,z,p ± ,Q 2 ) = M^ix.z.Q 2 ) , (6) 

as direct consequence of the kinematics of both terms in the evolution equations, 
eq.JSJl. To be definite we must say that the soft and collinear factorization of these 
objects is at the moment only a conjecture. 
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3 Solutions to evolution equations and phenomenology 



In our opinion the study of target fragmentation phenomenology is indeed the 
most intruiguing aspects offered by the novel equations we have proposed. However, 
due to lack of the aforementioned factorization theorem for TMD fracture functions, 
we think that current fragmentation phenomenology is a good benchmark for our 
evolution equations since a great amount of data to confront with is available. In 
what follows we thus describe the method for solving the evolution equations in 
the space- like case only, eq.©, since similar arguments applie also in the time- like 
case, eq. Solutions in both kinematics allows us to compare directly with current 
fragmentation phenomenology. 

As in the longitudinal case, distributions at a scale Q 2 > Qq are calculable if we 
provide a non-perturbative input density at some arbitrary scale Qq. We assume as 
initial condition a longitudinal distribution functions |12II13| times a cc-independent 
gaussian transverse factor: 

k 2 /<k 2 > 

T P (x B ,Q 2 ,k ± ) = F P (x B ,Q 2 ) C - . (7) 

7T < fej_ i > 

The gaussian fc^-distribution is used to model partons intrinsic momenta inside 
hadrons at the initial scale. This assumption is reasonable in the valence region 
x 6 0(1) but questionable in the gluon-dominated small-x range, a relevant issue 
for HERA kinematics, where the transverse factor, we argue, should become in 
some way x-dependent and also modified in the shape. Average transverse momen- 
tum is then set to < kj_ iq . >~ 0.25 GeV 2 . We choose a finite difference method 
for solving the evolution equations since it give us a direct control on the starting 
distributions and we set the splitting functions P&'(x) at lowest order (n — 0) in 
their expansion. The r.h.s. of the evolution equations are thus integrated on a bidi- 
mensional grid in (x,k±). At each iteration on the Q 2 variable, the normalization 
conditions, eq.(0J, is checked to reproduce ordinary longitudinal distributions. We 
solve the (2n/+l)— dimensional matrix equations in the space of quarks, antiquarks 
and gluons. The resulting space- like TMD distributions show a x— dependent broad- 
ening of the transverse momentum spectrum with increasing Q 2 . Furthermore the 
factorized form of eq.JTJ is not preserved under evolution. This result is expected 
since the transverse part of the TMD distributions, eq.©, mixes longitudinal and 
transverse variables. The spectrum tail shows a inverse-powerlike behaviour of the 
type (fcj_) _Q , with 1 < a < 2 depending on x, at high k\. Noticeable deviations 
from a guassian behaviour at incresing Q 2 are observed. The more pronunciated 
broadening in the transverse spectrum is more appreciable in the TMD gluon dis- 
tribution at small x. Such an effect is expected since P gq (x), P gg (x) ~ 1/x as x — > 0. 

As previosly stated, with TMD distributions at hand, we can directly recon- 
struct the structure functions .4(2:5, Zh,Q 2 , Ph±) |15) . in the current frag- 
mentation region according to 
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H 2 (xB,z h ,P hX ,Q 2 ) = e l J d 2 k ± d 2 P± 6 (2) (z h k ± +p ± - P h± ) ■ 

*=?> 9 

■ r P {x B ,A, k±, ) l£(z h> n%,p ± ) C(Q 2 , fi 2 F , fi 2 D ) , (8) 

where the factorization of eq.© in terms of TMD distributions has been explictely 
proven in ^1] and is accurate up to powers in (P%±/Q 2 ) n for transverse momenta 
\Ph±\ — Aqcu- The standard SIDIS variables are defined as xb — Q 2 /(2P ■ q) 
and Zh — [P ■ Ph)/{P ■ q) ■ Since we are using LL evolution equations, we set 
C = 1 and the factorization scales p? F — p? D — Q 2 . In our opinion is interesting 
to note that gluon, not directly coupled to the lepton current at LL in eq.lJSJ, 
partecipate the process through the evolution equations, eq.l© and eq.©. We are 
thus, even at lowest order, indirectly sensitive to the transverse spectra of the gluon 
which manifests its dynamical effects mostly on sea-quarks transverse spectra at 
small x. In the opposite kinematic region, x — ► 1, we are instead sensitive to soft 
gluon resummation effects which can effectively taken into account in our evolution 
equations setting a s (Q 2 ) — > a s ((l — x)Q 2 ) JHI- Once resummation effects are taken 
into account, the valence region it is the best suited kinematical region in which 
extract the non-perturbative quark intrinsic transverse momentum. 

Eq. (jHl can be shown to include also target fragmentated hadrons using TMD 
fracture functions In this case, as previously stated, a factorization theorem 
for TMD fracture functions is lacking and the factorization is only conjecturized. 

4 Conclusions 

We have shown that radiative transverse momentum can be quantitatively taken 
into account using TMD DGLAP evolution equations and that the resulting TMD 
distributions show a k\— smearing due to pQCD dynamics. The extension can be 
performed both in the current and in the target fragmentation region. Delicate 
issue concerning the small-x gluon distribution, impact of soft gluon resummation 
to gauge properly intrinsic transverse momentum of partons and peculiar feature 
of the transverse spectrum of TMD distributions are presented. The formalism we 
have developed can be extended to the polarized case both in current and in the 
target region [T7). Furthemore, if factorization holds for TMD fracture functions, 
we will have a multi-dimensional description of SIDIS cross-sections, valid onto the 
whole phase space of the produced hadron. A detailed derivation of our results and 
a complete phenomenological study will appear in a forthcoming paper |18) . 
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